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Considering the mixture after muon-catalyzed fusion (uCF) reaction as overdense plasma, we analyze muon 
motion in the plasma induced by a linearly polarized two-colour laser, particularly, the effect of laser parameters 
on the muon momentum and trajectory. The results show that muon drift along the propagation of laser and os- 
cillation perpendicular to the propagation remain after the end of the laser pulse. Under appropriate parameters, 
muon can go from the skin layer into field-free matter in a time period of much less than the pulse duration. The 
electric-field strength ratio or frequency ratio of the fundamental to the harmonic has more influence on muon 
oscillation. The laser affects little on other particles in the plasma. Hence, in theory, this work can avoid muon 
sticking to a effectively and reduce muon-loss probability in wCF. 
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I. INTRODUCTION 


The muon catalyzed fusion (uCF) of hydrogen isotopes, 
especially D-T fusion, has been studied as a realizable can- 
didate of an energy source as thermal energies [1—7]. In the 
liquid mixture of Dz and T2, the muon plays the role of a cat- 
alyzer through the formation of a muonic molecule. After the 
fusion process, the muon is released normally and again it is 
utilized for a subsequent fusion until the muon decays. How- 
ever, in fact, after the uCF, some muon may be strapped by 
a particles and form into (ay:)*. This has become a major 
limitation on the number of fusion cycles catalyzed by each 
muon [8, 9]. Because of muon-loss, CF is far from the aim 
to attain energy plus. In the past decade, researchers found 
sticking phenomena and new insights in jsCF [10-14], and 
in experiments wCF numbers per muon extend to 150 during 
the lifetime of ~2.2 x 1079s [15]. Althoght this catalyzed 
efficiency is not enough for atomic energy application, muon 
regeneration from the a-sticking has been increasingly stud- 
ied. 

In this work, we were trying to avoid muon-a sticking by 
inducting an external force. Here, the laser is a good choice. 
It is well-known that ultra-short, ultra-intense laser pulses 
are potentially useful for a variety of applications. Recently, 
ultra-short lasers have been applied extensively to studies of 
electron acceleration schemes due to their potential of offer- 
ing compact and low-cost setups through accelerator gradi- 
ents several orders of magnitude higher than those of con- 
ventional RF linear accelerators [16-18]. Unlike one-color 
laser, a two-color laser beam is the sum of two copropagat- 
ing pulsed radially polarized laser beams, with central angu- 
lar frequencies w and nw, of unequal pulse duration and peak 
power. Also, the relative phase of the two-color laser pulses 


* Supported by Natural Science Foundation of Scientific and Tech- 
nological Innovation Programs of Higher Education Institutions in 
Shanxi (No.2009 1025) 

* Corresponding author, shi-chunhua@ 163.com 


adds one more degree of freedom [19]. Previous simulations 
have shown that a one-color pulsed radially polarized laser 
beam can accelerate an initially stationary electron only up to 
40% of the theoretical energy gain limit [20, 21], while a two- 
color pulsed beam can accelerate an electron by over 9096 of 
the one-color beam's theoretical gain limit, for a given total 
energy and pulse duration [22]. So, in order to make the muon 
well out of the mixture, the choice of the two-color laser be- 
comes especially important [23, 24]. 

In this paper, the two-color laser is induced in the D-D fu- 
sion, and muon can be moved at a certain track so that pu is 
not captured by o. To study the muon drift of D-D fusion 
in the two-color laser, it is necessary to choose some appro- 
priate laser parameters. Here, we explore the possibility of 
the muon drift and trajectory, particularly; study influences of 
various laser-pulse parameters (i.e., frequency, intensity, and 
relative phase) on the muon drift and investigate coherent ef- 
fects on other particles in a two-level field. 


II. COMPUTATIONAL MODEL AND METHOD 


When a cycle of uCF(d-d): u+d+d — ddu + o--n- pis 
just finished, the a and p particles are existing independently, 
ie. the jj ions are not captured by the o particles. Then, at 
the moment, there is the mixture which consists of o, neu- 
tron (n), electron (e), d, and yz in the fusion reaction target. 
It is known that the use of a plasma medium is an attractive 
way of achieving laser-driven electron acceleration [17]. So, 
we consider the above mixture as the plasma with a frequency 
Wp, (We = we, + Wa. + Wu Wpas Wpe and Wp, are frequency 
of a, e and p, respectively). Also, as wpe > Wpu > Wpa, We 
judge w? = we. 

At the same time of the above mixture formation, the two- 
color laser pulse is turned on and off adiabatically, i.e. 


E(t) = f(t)- {Erexp(—iwet) + En exp[—i(wnt + v)]), (1) 


where, subscript f denotes the fundamental part in the two- 
color laser pulse and subscript h denotes the harmonic part; 
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E; and Ey are amplitude of the electric-field strength of the 
fundamental laser field and the harmonic field, respectively; 
for the frequency wf and wp, the n' harmonic field means 
Wh = nwr; is the relative phase between the fundamen- 
tal field and the harmonic; and f(t) is the shape of the laser 
pulse, i.e. 


f(t) = sech?(t/r), (2) 


where, 7 is the pulse duration. When the plasma frequency 
Wp meets the relation of wp 2» wn > w (w = (wr + wy) /2 is 
the equivalent frequency of a two-color laser), the laser beam 
can be considered propagating in overdense plasma [25]. 

For the sake of simplicity, we assume that the linearly po- 
larized laser beam propagates perpendicular to the plane sur- 
face of the target. Let us direct the axis X along the prop- 
agation of the laser pulse, the axis Y along the electric-field 
strength, and the axis Z along the magnetic-field strength. 
Throughout the paper, the system of units m, = e = c = 1 
is used. A laser beam propagates in overdense plasma with a 
plasma frequency wp greater than the laser frequency w in the 
skin layer having a depth ô = 1/wp. The boundary conditions 
for the electric and magnetic fields on the surface vacuum- 
overdense plasma are of the well-known form: 


in 2 2 
p"- T D RI ver 

im. _2VE a 
LEE Rew e 2F(t), 


where, F' is the electric-field amplitude inside the plasma, 
F(t) meets the relation of F(t) = |E(t)| and varepsilon is 
the dielectric constant of plasma, ie. € = 1 — ws fur? fd 
—w? fu. 

Then, the Newton classical equations for muon momentum 
P, and P, inside the skin layer (x > 0) are of the form: 


pers =2v, f (t)exp(—a/d){ E sin(wt) 


d 
+ Ey sin(w,t + y)} 
kd =2— f(texp(—x/3){Eycos(wyt) (a) 
Wp 


+ En cos(wpt + y)} 
— 2v, f (t)exp(—2/d){ Ey sin(wt) 
+ Ej sin(wypt + p)} 


where k = m, /m. e 207, 6 = 1/wp is the depth of the skin 
layer, and x = 0 is the surface of the target. The velocity v; 
and vy are given by 


da P; 
Vy = = ; 
dt 12 + P2 + P2 
T y 
dy P, (5) 
V, = — 
” dt 12 4 P24 P2 
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Ill. NUMERICAL RESULTS AND ANALYSIS 


We assume that there is no muon motion along the 
magnetic-field strength, i.e. the axis Z. Besides, collisions 
of muons are neglected in our consideration. Eqs. (4) and (5) 
can be solved with the initial conditions that at the time in- 
stance t = —oo a muon is at rest when x = y = 0. Of course, 
the average muon velocity along the Y axis is zero. After 
the end of a laser pulse, a muon has zero velocity in the Y 
direction and nonzero velocity in the X direction. 


A. The effect of electric-field strength 

The relativistic plasma frequency is wp = (4mn./y)'/?, 
y = [1 + (F9 /w)?]"/? is the relativistic factor. From Eq. (1) 
to Eq. (5), we fixed the value w, /w, = 21/2 and wet /m = 40 
in our numerical computation which correspond to the pulse 
duration 7 = 83fs and the electron number density ne = 
8 x 107? cem ?. For the two-color laser, we take the 2™ har- 
monic field wp = 2w, and the relative phase y = 0. Based on 
the above values, by solving the Eqs. (4) and (5), the muon 
motion is determined only by one dimensionless parameter 
g = 4Er/(3wr) and the amplitude ratio of the electric-field 
strength Er/ E, . The value g = 1 corresponds to the peak 
laser intensity of 5 x 10!" W/cm? (at the laser photon energy 
hw = 1.5eV). 
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Fig. 1. The dependences of the muon momentum P; on the dimen- 
sionless parameter g. 


Figure | shows the dependence of muon momentum P, (in 
units of mec, in this paper, the units of the momentum are all 
mec.) on the value of g and the different ratio of Er to Enp, 
where we take E, = E; and 3E, = Ep, respectively. The 
value g from 0 to 100 corresponds to the peak laser inten- 
sity of 0 ~ 5 x 10?! W/cm?. The results show that the two 
curves change in a similar way basically. The momentum 
increases with the laser intensity, but it is a nonrelativistic 
quantity (P, < mec) even at the relativistic peak laser in- 
tensity 5 x 10?! W/cm?. Fig. 1 also shows the different ratio 
of the fundamental electric-field strength to the harmonic has 
effects on the muon momentum P, but the effect is obvious 
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Fig. 2. (Color online) The muon momentum P; and P, change for the value of the dimensionless parameter g = 80. 


at g = 0 ~ 70 and is not obvious at g > 70. 

Figure 2 demonstrates the muon momentum P, and P, at 
the leading edge of the laser pulse for the typical value of 
g = 80 (the peak laser intensity of the fundamental wave 
is 3.2 x 10?! W/cm? according to Eqs. (4) and (5), with 
Ey, = Erin Fig. 2(a1-a3) and 2E, = Er in Fig. 2(b1-ba). As 
can be seen in the figures, a muon goes from the skin layer 
inside field-free matter before the laser pulse (ox sec h?(t/7)) 
reaches its maximum because of the small depth of the skin 
layer. i.e. a muon stops acceleration on the rear border of the 
skin layer. It is also obvious that P,—0 when a muon leaves 
the skin layer, because of adiabaticity of the laser pulse. In ad- 
dition, the different E;/ E ratio affects little on the muon mo- 
mentum in the X axis, but obviously in the Y axis (Figs. 2(a3) 
and 2(b3)). 

The muon trajectory in the plane for g = 80 is shown in 
Fig. 3. The dimensionless muon coordinates x and y are 
given in units of the skin depth 6. The amplitude of muon 
oscillations in the transverse direction Y is much less than its 
drift motion in the direction of laser pulse propagation (the .X 
axis). This indicates that a muon penetrates many skin depths 
into field-free matter even before the laser pulse reaches its 
maximum. The amplitude and shape of muon oscillations 
in the transverse direction Y in Figs. 3(a) and 3(b) differ 
from each other, due to different harmonic wave electric-field 
strength. At a greater m = Ert/ Ey, the amplitude of muon os- 
cillations in the transverse direction Y reduces and the shape 
is more complicated. 

To explain that a muon goes from the skin layer inside 
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Fig. 4. (Color online) The muon motion in one-color laser field. 


field-free matter or penetrates many skin depths into field- 
free matter before the laser pulse reaches its maximum, we 
make an estimate qualitatively. We consider that the muon 
drift velocity in the direction of laser pulse propagation has 
the estimate v, = N(Ert/wy)? from Eq. (4). Then, the time 
tp for penetration of the skin layer having the depth ô is ap- 
proximately, ie. tp = 5/vz = wy/(NE2) ~ 48fs. The 
result shows £j < 7, ie. the time tp at which a muon goes 
into field-free matter from the skin layer is much less than the 
laser duration 7. Thus, a muon goes quickly from the skin 
layer into field-free matter. 


Besides, in order to better illustrate effects of the two-color 
laser intensity on muon motion, we calculate the muon mo- 
mentum and trajectory in one-color laser (Fig. 4, in a peak 
laser intensity of 6.4 x 10?! W/cm? and other parameters be- 
ing the same value as the above. In Fig. 4(c) the maximum 
of muon momentum P, in the one-color laser field is less 
than in Fig. 2(a,). It shows that the influence of the two-color 
laser intensity on the momentum is not a simple addition of 
two one-color laser beams, and a two-color pulsed beam can 
accelerate muon to the same maximum with lower intensity. 
In Figs. 4(a) and 4(b), the oscillations of muon momentum 
P, in the one-color laser field are different from Fig. 2(a2) 
and 2(a3). The muon oscillations have just one peak in ev- 


ery cycle in one-color laser field (Fig. 4(b)) but two peaks 
in two-color laser field (Fig. 2(a3)). The muon trajectory in 
the one-color laser field is shown in Fig. 4(d), where the am- 
plitude of muon oscillations in the transverse direction Y is 
obviously less than that in the two-color laser field (Fig. 3(a)). 
In Fig. 3(a), the muon oscillation in each cycle has two peaks 
(red and black arrows) in the two-color laser field, but in 
Fig. 4(d) it has a single peak (red arrow). Certainly, in the 
one-color laser field, a muon goes from the skin layer inside 
field-free matter before the laser pulse reaches its maximum, 
and P,—0 when a muon leaves the skin layer. However, a 
muon takes a longer time to go into field-free matter from the 
skin layer in the one-color field than in the two-color field. 


B. The effect of frequency 


To evaluate the frequency effect, we still fix the pulse dura- 
tion 7 = 83 fs and the electron density n, = 8 x 107? cm^?. 
For the two-color laser, we take 2E, = Ey, Eg = 
8.0 x 102° W/cm? and o = 0. Therefore, the muon motion is 
determined only by the dimensionless parameter wp/w (w = 
(toe + wp) /2) and the frequency ratio of the fundamental laser 
field we to the harmonic field wp. 
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Figure 5 shows the dependence of muon momentum Py 
on wp/w, and on wr/wp, at wn = We, Wn = 2wf and wn = 
Aur, respectively. The wp/w changes from 5 to 20 so that the 
plasma frequency wy is far larger than the laser frequency w. 
In Fig. 5, the three curves change in a similar way basically, 
but different frequency ratio affects the muon momentum P, 
especially at wp/w = 6 ~ 14. Besides, the three curves have 
different peak positions. 

Figure 6 shows the muon momentum P, and P, at the lead- 
ing edge of the laser pulse for typical value of the peak laser 
intensity of the fundamental wave being 8.0 x 10?? W/cm? 
and wp/w = 8 according to Fig. 5. We take wn = 2wr in 
Figs. 6(a1) and 6(a2) and wp = 4wp in Figs. 6(b1) and 6(b2). 
From the figures, a muon goes from the skin layer inside field- 
free matter before the laser pulse reaches its maximum, and 
P, — 0 when a muon leaves the skin layer. In addition, 
the frequency ratio affects the muon momentum obviously in 
the X axis and the Y axis. At a higher frequency ratio, i.e. a 
greater n = wy /ur, the maximum of the muon momentum P, 
increase, and so does the amplitude of the muon oscillations 
of P, in the Y axis, with a more complicated shape. 

Figure 7 shows the muon trajectory in the plane for the 
peak laser intensity of fundamental wave 8.0 x 107? W/cm? 
and wp/w = 8. The dimensionless muon coordinates x and y 
are given in units of the skin depth ó. The amplitude of muon 
oscillations in the transverse direction Y is much less than its 
drift motion in the direction of laser pulse propagation (the X 
axis). This indicates that a muon penetrates the skin depths 
into field-free matter even before the laser pulse reaches its 
maximum. At a greater n = w/w, with an increased muon 
momentum in the .X and the Y axises, the amplitude of muon 
oscillations in the transverse direction Y increases while the 
frequency of muon oscillations reduces. Thus, the muon os- 
cillations in the Y axis reduce and at the same time the muon 
leaves the skin layer by a faster speed in the X axis. 

In the same way, based on above parameters, the time tp at 
which a muon goes into field-free matter from the skin layer 
is about 44fs. A smaller-than-7 tp means that a muon pen- 
etrates the skin depth into field-free matter before the laser 
pulse reaches its maximum. So, the muon leaves the skin 
layer faster to enter the free field, at wp/w = 8 and peak 
laser intensity of E; = 8.0 x 102° W/cm?, than the results in 
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Sec. III. A at much larger parameters of g and Fy. 


C. The effect of relative phase 


For the two-color laser with the parameters of wp/w = 8, 
E; = 8.0 x 107? W/cm?, electron number density n, = 
8.0 x 107? cm7?, 2E, = Er and wp = 2wp, the muon mo- 
tion is determined only by the pulse duration 7 and the rela- 
tive phase y. The relative phase in two-color lasers means the 
propagation distance in the plasma media. It can be converted 
into the media length L = A,/(27). 

Figure 8 shows the muon momentum P at the relative 
phase of y = 7/6 (L = àn/12) and the pulse duration 
T = 10 ~ 100fs. At 7 < 50fs, the muon momenturm P, is 
less than those in Sec. III. A and Sec. III. B, hence a longer tp, 
the time for penetrating the skin layer of the depth ó. Then, 
probably, once the laser pulse duration is less than 50 fs, the 
muon cannot penetrate the skin layer at the end of a pulse. 
Thus, to study the effect of the relative phase on the muon 
trajectory, we take 7 = 50 fs. 

Figure 9 shows the muon momentum P, at the leading 
edge of the laser pulse, with peak laser intensity of the fun- 
damental wave being 8.0 x 107° W/cm? at w,/w = 8,7 = 
50 fs, and the relative phase y = 0 (L = 0), 7/4 (L = An /8), 
71/2 (L = Xn/4) and 37/4 (L = 3A4/8). We also find that 
the relative phase y has obvious effects on the muon momen- 
tum in just the Y axis. The muon oscillations shape and the 
peak change with the relative phase. The red and the black 
arrows indicate peaks that change constantly with the rela- 
tive phase. Because the relative phase has less effect on the 
muon momentum P,,, the muon trajectories are similar to the 
results above. The muon speed in the X axis is not affected 
and the muon still can penetrate skin depths into field-free 
matter. However, as shown in Fig. 10, in the transverse direc- 
tion Y, the muon oscillations are in different amplitudes and 
shapes at different relative phases. In Figs. 9 and 10, at greater 
the relative phases, the muon oscillation peaks marked with 
red arrows increase while the muon oscillation peaks marked 
with black arrows reduce consistently. 


D. The effects of muon initial velocity 


For all above analysis, initial velocity of muon is neglected 
in our consideration, i.e. Eqs. (4) and (5) can be solved with 
the initial conditions that v(t = —oo) = 0. In fact, after every 
LCF reaction, velocity of muon is not zero but its value is not 
unique. Most muons are free, and their velocity depends on 
the energy release of nuclear fusion [13]. Some muons stick 
to alpha particles and form muonic helium ions, and the initial 
sticking probability is 0.91296 [26]. Where muonic helium 
ions are formed with energy of 3.47 MeV, they are slowed 
down toward thermal energy by collision with the surround- 
ing Də and DT molecules [27]. So, muonic helium ions can 
be stripped as a result of collisions and some muons regain 
free, probability of which is 25% [28]. Their velocities are 
much less. Other muonic helium ions may be stripped by the 
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Fig. 6. (Color online) The muon momentum P, and P, change for the value of the dimensionless parameter wp /w = 8. 
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Fig. 7. (Color online) The muon trajectory. 


30 40 50 60 70 80 90 100 
t (fs) 


10 20 


Fig. 8. The dependences of muon momentum P on the pulse dura- 
tion 7 at the relative phase of y = 7/6. 


external force [29] and velocity of free muon are related to the 
force. Based on the energy release of nuclear fusion, we take 
two value of initial velocity and solve the Eqs. (4) and (5). 


Figure 11 shows the muon momentum P, and P, at the 
leading edge of the laser pulse and where values of relevant 
parameters are the same with Fig. 2(a;—a3). The only differ- 
ent is the value of muon initial velocity. According to Eqs. (4) 
and (5), we take the initial velocity vs = vy = 0.001 in 
Fig. 11(a1-a3), and v, = vy = 0.0001 in Fig. 11(b1—b3), 
instead of v, = vy = 0 in Fig. 2(a1—a3), the unit of veloc- 
ity is c = 3 x 10° m/s. The results shows that a muon goes 
from the skin layer inside field-free matter before the laser 
pulse reaches its maximum, and P, — 0 when a muon leaves 
the skin layer. Comparing Figs. 11(a3) and 11(b3), the detail 
changes of P}, the muon initial velocity has obvious effect on 
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Fig. 10. (Color online) The muon trajectory. 


the muon momentum in the Y axis, while it has a little effect on the muon momentum in the X axis. In the case of super- 
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Fig. 11. The muon momentum P, and P, at g = 80 in muon initial velocity of 0.001c and 0.0001c (c = 3 x 10? m/s). 


intense laser pulses, the initial velocity in the X axis plays 
a little role in muon drift, for muon has relatively high drift 
velocities induced by the laser in the X axis. 


IV. DISCUSSION AND CONCLUSION 


In addition, with the induction of the two-color laser, we 
have to consider the influence of the laser on other particles 
in the overdense plasma. According to similar methods, we 
study the other particles trajectory in different muzzle veloc- 
ities. Because the muzzle velocity of other particles is small, 
the laser has little influence on them, i.e. the motion of those 
particles produced by the superintense laser pulse is far less 
than the motion of muons. Also, the duration of laser pulse 
is so short that the work does not introduce the disadvan- 
tage into wCF reaction in the superintense laser, but make the 
muon penetrate skin depths into field-free matter to avoid the 
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